09) 



J 



(12) 



(43) Date of publication: 

05.02.1997 Bulletin 1997/06 

(21) Application number: 95309127.9 

(22) DateoHiling: 14.12.1995 



Europfilsches Palentamt 
European Patent Office 

Office europeen des brevets (11) EP 0 757 243 A1 

EUROPEAN PATENT APPLICATION 

(51) int. CI 6 : Q01N 21/35 



(84) Designated Contracting Slates: 

AT BE CH DE DK ES FR GB GR IE IT LI LU MC NL 
PTSE 

Designated Extension States: 
LT LV SI 

(30) Priority: 31.07.1995 US 509696 
24.10.1995 US 509696 

(71) Applicant: INSTRUMENTATION METRICS, INC. 
Chanhassen, Maine 55317 (US) 



(72) Inventor: Schlager, Kenneth J. 
Elm Grove, Wisconsin 53122 (US) 

(74) Representative: Cross, Rupert Edward Blount et 
al 

BOULT WADE TENNANT 
27 Fumival Street 
London EC4 A 1PQ (GB) 



(54) Liquid correlation spectometry 

(57) A number ol optically active elements are pro- 
vided that enable the performance of accurate, depend- 
able chemical analysis in a medium having varying 
matrix backgrounds and substantial component interfer- 
ences. Optical transfer cells are provided that include 
positive correlation filter systems adapted to accept 
near-IR radiation and to pass a wavelength from that 
radiation having high correlation with a particular ana- 
lyte. A near-IR spectrometer apparatus for the measure- 
ment of anaiyte concentration in a complex medium is 
also provided. The apparatus features an optical trans- 
fer cell constructed according to the invention. The 
apparatus can be used in the non-invasive determina- 
tion of blood anaiyte concentrations in a variety of medi- 
ums. Methods of measuring the concentration of a 
blood anaiyte in a mammalian subject are also pro- 
vided. The methods entail the use of reflective near-IR 
radiation to interrogate a medium to obtain spectrally 
modified radiation. The radiation thus obtained is 
passed through an optical transfer cell constructed 
according to the invention to obtain a wavelength having 
high correlation with an anaiyte of interest and substan- 
tially no correlation with interfering constituents present 
in the medium. 
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Description 

r.rn ^.Ra f ? mnr.B to Re lated Applications 

5 This application is a continuation-in-part of U.S. Patent Application Serial No. 08/509,696, filed July 31,1 995, from 

which priority is claimed pursuant to 35 U.&C. §120, and which disclosure is hereby incorporated by reference in its 
entirety. 

Technical Field 

10 

The present invention relates to a spectrometer device for quantifying target chemical analytes in a complex aque- 
ous medium and methods ot use thereof. The invention finds application in a wide range of chemical analyses, and par- 
ticularly in non-invasive spectrophotometry analysis of Wood analytes. One important application of the invention 
involves the precise measurement of Wood glucose using non-invasive sampling techniques. 

is 

Barkynund of the Invention 

As patient out-care services become more prevalent in modern health care, the need for low-cost reliable instru- 
mentation becomes more acute, particularly for applications where blood analyte chemistry monitoring must be carried 
20 out. One such application is in home-monitoc ing of Wood glucose levels by diabetics. 

Diabetes is a major health concern in the United States where it has been estimated that there are upwards of 13 
million persons afflicted with the condition. Lipsett, L (1993) "Statistics; Prevalence, Incidence, Risk Factors and Com- 
plications of Diabetes," American Diabetes Association Bulletin April 9, 1993. Treatment of the more severe form ol the 
disease, Type I insulin-dependent diabetes (where the pancreas secretes no insulin) requires one or more insulin injec- 
25 tions per day. Insulin controls the utilization of glucose or sugar in the Wood and prevents hyperglycemia which, if left 
uncorrected, can lead to ketosis, coma and death. 

Hyperglycemia in diabetics has been correlated with several long-term effects of diabetes, such as heart disease, 
atherosclerosis, blindness, stroke, hypertension and kidney failure. Normally, a persop's blood glucose level may vary 
between 60 and 100 milligrams per deciliter, whereas in a diabetic, the levels may vary from 40 to 50X5 milligrams per 
30 deciliter. Such drastic blood glucose level swings must be avoided in order to prevent the serious complications of the 
disease. Accordingly, blood glucose concentrations need to be periodically monitored, in some instances up to four 
times daily, so that diabetics may monitor their Wood glucose level and vary caloric intake, diet and insulin accordingly. 

Traditional blood glucose monitoring methods require the diabetic to draw a blood sample (e.g., by a finger-tip 
lance) for each test and to read the glucose level using a glucometer (a spectrophotometer that reads glucose ooncen- 
35 tralions) or a colorimetric calibration method. Such invasive blood extractions create a painful and tedious burden to the 
diabetic and expose the diabetic to the possibility of infection, particularly in light of the frequency of testing which is 
necessary. These considerations can lead to an abatement of the monitoring process by the diabetic. 

The value of frequent monitoring of blood glucose as a means to avoid or at least minimize the complications of 
Type I diabetes has been well established. Further, patients with Type II (non-insulin-dependent) diabetes can benefit 
40 from blood glucose monitoring in the control of their disease, however, the inconvenience and traumatic nature of 
known procedures discourages such monitoring in many patients. 

Accordingly, there is a recognized need in the art for a simple and accurate method and apparatus for non-inva- 
sively measuring blood analyte concentration, particularly in the context of blood glucose monitoring by diabetics. Var- 
ious types of apparatus and related methods for non-invasive measurement of glucose concentration and other blood 
45 analyte concentrations have been described in the art Such devices use traditional methods of near infrared (near-IR) 
analysis based on the measurement of absorbance at one or more specific wavelengths. 

U.S. Patent No. 5.267,152 to \ang et al. descrfoes a non-invasive technique of measuring Wood glucose concen- 
tration using near-IR radiation diffuse-reflection laser spectroscopy. The device uses only a narrow portion of the IR 
spectrum referred to as the "water transmission window" (those wavelengths between 1300 and 1900 nm). The optl- 
so caJly controlled light is impinged on a tissue source and then collected by an integrating sphere. The collected light is 
then analyzed and blood glucose concentration is calculated using a stored reference calibration curve. 

US. Patent No. 5,086,229 to Rosenthal et al. describes a device for measurement of blood glucose using analysis 
of near-IR energy following interaction with venous or arterial blood, or following transmission through a blood contain- 
ing body part. A transmissive technique is described wherein at least one source of near-IR energy (of a selected band- 
55 width) is provided for passage through a sample. Measurement of glucose concentration is carried out using 
quantitative analysis algorithms which generally employ ratios and gaps to create mass coefficients of absorption. 

U.S. Patent No. 4,975,581 to Robinson et al. describes a method and apparatus for determining analyte concen- 
tration in a biological sample based on a comparison ot infrared energy absorption (i.e., differences in absorption at 
several wavelengths) between a known analyte concentration and a sample. The comparison is performed using partial 



2 



4 



EP0 757 243 A1 

least squares analysis or other multivariate techniques. 

U.S. Patent No. 4,882.492 to Schlager describes a method and apparatus for non-invasive determination of Wood 
analyte concentrations. Modulated IR radiation is directed against a tissue sample (e.g., an ear lobe) and either passed 
through the tissue or impinged on a skin surface where it is spectrally modilied by a targel analyte (glucose). The spec- 

5 trally modified radiation is then split, wherein one portion is directed through a correlation cell and another through a 
reference cell. Intensity of the raolation passing through the cells are compared to determine analyte concentration in 
the sample. 

Additionally, similar devices have been described for use in determination of analyte concentrations in complex 
samples. 

io U.S. Patent No. 5,252,829 to Nygaard et al. describes a method and apparatus for measuring the concentration ot 
urea in a milk sample using an inlrared attenuation measuring technique, wherein milk samples (of substantially known 
and constant component concentration) are analyzed for IR attenuation at discrete wavebands where urea, fat lactose 
and protein absorb, Multivariate techniques are carried out to determine spectral contributions of each known compo- 
nent using partial least squares algorithm, principal component regression, multiple linear regression or artificial neural 

75 network learning. Calibration is carried out by accounting for the component contritions that block the analyte signal 
of interest (urea). Thus, Nygaard et al. describe a technique of measuring multiple analyte infrared attenuarjons-and 
then compensating for the influence of background analytes on the measurement ot a selected analyte (e.g., compen- 
sating tor background contribution to the urea signal). 

U.S. Patent No. 5,242,602 to Richardson et al. describes methods far analyzing aqueous systems to delect multiple 

20 active or inactive water treating components. The methods involve determination of the absorbance or emission spec- 
trum of the components over the range of 200 to 2500 nm. and application of chemometrics algorithms to extract seg- 
ments of the spectral data obtained to quantify multiple performance indicators. 

U.S. Patent No. 4,306,152 to Ross et al. describes an optical fluid analyzer designed to minimize the effect of back- 
ground absorption (the overall or base level optical absorption of the fluid sample) on the accuracy of measuring a com- 

25 ponent in a turbid or otherwise hard to analyze liquid sample. The apparatus measures an optical signal at the 
characteristic optical absorption (of a sample component of interest) and another at a wavelength selected to approxi- 
mate background absorption and then subtracts to reduce the background component of the signal. 

None of the aforementioned devices or methods have been capable of providing low cost, accurate determination 
of Wood analyte concentration in field or home-testing applications. Particularly, known devices have not provided tor 

3d the accurate determination of Wood analyte concentrations in the midst of multiple interferences from other analytes 
that are active in the same spectral region as the analyte of interest. 

SMrnmaryolthBlriYgntiQn 

6 Accordingly, it is a primary object ot the invention to address the above-mentioned needs in the art, by providing 
optically active elements that enable one to perform accurate, dependable chemical analysis in an aqueous medium 
with varying matrix background and substantial component interferences. The optically active elements allow for opti- 
mal analyte measurement at low cost tor a wide variety of medical, industrial and environmental applications. 

Spectrometers are generally known for use in a wide variety of areas such as medidna scientific research, indus- 

40 trial applications and the like. In one aspect of the invention, an optical transfer cell is provided for use in spectrometer 
devices which can be used to quantify target chemical analytes in an aqueous medium. In this regard, the optical trans- 
fer cell is particularly intended for use in measurements where there exist complex spectral backgrounds (e.g., where 
interfering substances introduce significant measurement error), and especially where concentrations of background 
substances greatly exceed those of a target analyte. The subject optical transfer cell finds application in a wide range 

45 of chemical analyte analyses carried out in a diversity of aqueous media. 

In one embodiment, an optical transfer cell is provided which comprises a positive correlation filter system. The 
subject filter system is adapted to accept near-IR radiation and to pass a wavelength from that radiation which is 
selected to have a high correlation with a particular analyte. The optical transfer cell can thus be used in the measure- 
ment of analyte concentrations In complex aqueous media, such as in the measurement of Wood analytes in body tis- 

50 sue. In another, related embodiment, a positive correlation filter system is provided that is capable of accepting near-IR 
radiation and selectively passing a plurality of wavelengths having high correlation with a particular analyte. 

In an alternative embodiment, a positive correlation filter system is provided that has been formed from a plurality 
of filter layers, each layer having a selected optical activity (absorption characteristic), whereby the system is capable 
of passing a population of wavelengths selected to have high correlation with a particular analyte concentration. In one 

55 particularly preferred embodiment, the aforementioned positive correlation filter system has at least one layer that 
includes wavelength weighting factor means. In this manner, the wavelength weighting factor provides enhanced corre- 
lation between a passed wavelength and analyte concentration in a particular medium. 

Thus, in order to overcome absorption interference caused by background components in a given sample, the 
present invention provides a number of positive correlation filter systems which are designed to provide a high correla- 
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tion with a target analyte and significantly tower correlations with other background absorbing materials. That is. posi- 
tive correlation titration relates to the transmission of a pattern of light through an optically active lilter which selectively 
passes "weighted" wavelengths; wherein the passed wavelengths correlate closely with a target analyte concentration 
in a given sample. 

5 In another aspect of the invention, positive correlation filter systems are synthesized using mathematical optimiza- 

tion techniques such as least squares analysis, linear rotated principal component analysis (RPC), or nonlinear genetic 
algorithm analysis. Thus, weighting factors can be determined by application of the aJbove-descrked optimization tech- 
niques to an original absorption spectrum of a selected analyte. In preferred embodiments, positive correlation filters 
are formed from a plurality of layets. each layer having an appropriate absorption characteristic supplied by a weighting 

to factor, wherein the filter selectively passes a population of wavelengths are weighted according to positive correlation 
to target analyte concentration. 

In yet a further aspect of the Invention, an apparatus for the accurate measurement of analyte concentration In 
complex aqueous media is provided. The apparatus includes a source of near-IR radiation, a means for reflectively 
transmitting near-IR radiation from the source to an aqueous medium which is to be interrogated, and a means for 

is retransmitting spectrally modified radiation reflected from the aqueous medium. The apparatus also includes means for 
receiving and directing the spectrally modified radiation into a beam path which is communicated to an optical transfer 
cell disposed within the beam path. The transfer cell comprises a positive correlation filter constructed according to the 
invention. The subject positive correlation filter is adapted to selectively pass a wavelength from the beam of spectrally 
modified radiation, wherein the wavelength has high correlation with a selected analyte concentration. The passed 

20 wavelength is then communicated to a means for receiving said information and converting the same into a signal which 
is representative of the intensity of the wavelength (e.g.. a photodetector or a near-IR photodiode array), and that signal 
is transmitted to a means for calculating analyte concentration. 

In the practice of the invention, the above-described apparatus can include any of the positive correlation filter sys- 
tems as described herein. In particular embodiments, the subject devices can be used to measure blood analyte con- 

25 centration in body tissue. The blood analyte to be measured can be selected from the group of calcium, potassium, 
sodium, chloride, bicarbonate (C0 2 ), hydrogen ion concentration (pH), glucose, urea (blood urea nitrogen, or BUN), 
hematocrit and hemoglobin, The source of near-IR radiation preferably supplies radiation including at least wavelengths 
in the range of approximately 1 100 to 2500 nm. In one particularly preferred embodiment, the apparatus can be used 
in the non-invasive analysis ol a blood analyte such as blood glucosa 

30 In an additional aspect of the invention, a non-invasive method of measuring the concentration of a blood analyte 
in a mammalian subject is provided. The method includes the steps of: 

(a) directing near infrared radiation to a source containing the blood analyte; 

(b) causing the radiation to be reflected by the source to provide a spectrally modified beam; 

3S (c) directing the modified beam into a path through an optical transfer cell disposed in the beam path. The optical 
transfer cell includes a positive correlation filter means that is adapted to accept the spectrally modified beam and 
to pass a wavelength from that beam, wherein that wavelength has been selected to have high correlation with the 
analyte concentration: 

(d) measuring the intensity of the passed wavelength and generating a signal representative of the intensity of the 
40 wavelength: and 

(e) converting the signal into an indication of analyte concentration. 

In one case, a non-invasive technique for measuring a Wood analyte is provided which enables reflective, transcu- 
taneous tissue measurement of Wood glucose in the forearm area. Such measurement gives rise to a more uniform 
4$ representation of blood glucose than previously possible using transmissive measurements in the finger area. 

Brief Description of the Figures 

Figure 1 is a diagrammatic representation of one embodiment of a correlation spectrometer constructed according 
so to the invention. 

Figure 2 is a diagrammatic representation of an alternative embodiment ol a correlation spectrometer constructed 
according to the invention. 

Detailed Description of the Preferred Embocfiments 

55 

Before the invention is described in detail, it is to be understood that this invention is not limited to the particular 
component parts of the devices or methods described, as such may vary. K is also to be understood that the terminol- 
ogy used herein is tor purposes of descrbing particular embodiments only, and is not intended to be limiting. It must be 
noted that, as used in the specification and the appended claims, the singular forms "a," "an" and "the" include plural 
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referents unless the oontext dearly dictates otherwise. Thus, for example, reference to "an analyte" includes mixtures 
of analytes, reference to "an optical transfer cell" includes two or more such optical transfer cells, "a means for reflec- 
tively transmitting radiation" includes two or more such means, "a wavelength" includes two or more wavelengths, "a 
chemometrics algorithm" includes two or more such algorithms, and the like. 
5 In this specification and in the claims which follow, reference will be made to a number of terms which shall be 

defined to have the following meanings: 

As used herein, a "reference sample" refers to an aqueous composition comprising a known concentration of an 
analyte of interest. A "standard" or "calibranf is a reference sample used for establishing the response of a measure- 
ment system to an analyte. As such, a calibration reference is typically any of the standards of various types that indi- 
ro cate whether an analytical instrument is working within prescribed limits and by which adjustments may be made to the 
analytical measurement system to correct for any deviation from the prescribed limits. 

"Calibration" refers to the process of relating an instrumental response to a concentration or property of a targeted 
matrix with the objective of developing a means to provide accurate and precise predictions in the future. The calibration 
step produces a direct interpretation of a measurement through a mathematical transformation that models the relation- 
is ship between a concentration and a spectral response. 

Methods of calibration are divided into estimation and classification techniques. "Estimation techniques" relate to 
the process of predicting the level of analyte concentration, while "classification techniques" determine which one of 
several predefined groups a sample will belong to. Methods of estimation are particularly useful in applications wherein 
the measurement of concentration is desired. To calibrate an instrument for estimation, it is necessary to have a set of 
20 exemplary measurements with known concentrations referred to as the calforation set (e.g., reference set). This Get 
consists ot m samples, each with n instrument variables contained in an m by n matrix (X), and an m by 1 vector (y). 
containing the concentration. If a priori information indicates the relationship between the measurement and concen- 
tration is linear, the process of calibration will attempt to determine an n by 1 transformation or mapping (£>), such that 

25 y 3 *b 

is an optimal estimate of y accorcfing to a predefined criteria. Estimation techniques are known in the art and include 
multiple regression, partial least squares regression, principal component regression and rotated principal components 
regression. 

30 The classification of a measurement is the primary focus of pattern recognition. Assuming a calibration set X. in 
which the class of each sample is known, methods of classification are used to develop a calibration model that will 
allow the classification of unknown measurements. Methods of classification are generally known in the art. and include 
template matching which involves the use of templates and a measurement of similarity, and discriminant analysis that 
optimizes the separation of samples in a given class through a rotation of measurements. 

35 The general validity of a calibration model can be verified by estimating the concentrations of independent meas- 
urements, a "test setr that has not been used in the calibration procedure. When the number of samples available for 
calibration are limited, the data set can be segmented, wherein part of the data set is used as the training set, and the 
remaining data used for the test set. 

"Chemometrics" relates to the application of mathematical, statistical and pattern recognition techniques in chem- 

40 ical analysis applications. See, e.g., Brown et al. (1990) Anal. Chem. 62:84-101 . Chemometrics is practiced herein in 
the context of developing and using non-invasive medical diagnostic instrumentation that employs advanced signal 
processing and calibration techniques. Signal processing is used to improve the accessibility of physically significant 
information in analytical signals. Examples of signal processing techniques include Fourier transformation, first and 
second derivatives, and digital or adaptive filtering. In the context of chemometrics, "calibration" refers to the process of 

45 relating data measurements to a chemical concentration for the purpose of quantification. Particularly, statistical cali- 
brations using chemometric methods can be used to extract specific information from a complex set of data. Such meth- 
ods of calibration include linear regression, multiple-linear regression, partial linear regression, and principal 
components analysis. In other applications, calibrations can be carried out using artificial neural networks, genetic algo- 
rithms and rotated principal components. 

so Instrumentation that detects information for one or more constituents in a complex chemical matrix must rely upon 
analysis algorithms (such as those derived using chemometrics) in order to reveal information thai is specific for one or 
more chemical constituent. Chemometric techniques can be used to compare unknowns with calibrated standards and 
data bases to provide advanced forms of cluster analysis, and to extract features from an unknown sample that can be 
used as information in statistical and mathematical models. 

ss In the practice of the invention, quantitative determinations of chemical analyte concentration using information 
from UV-vis-NIR spectroscopy can be made using computer algorithms, the parameters of which are determined using 
chemometric calibration sequences termed "learning sets." Learning sets are generally based on a large number of 
known (e.g., reference) samples that are used to determine the parameters of the algorithm. The number of samples 
that are required in any one learning set will be a function of the complexity of the sample matrix and of the number of 



5 



EP 0 757 243 A1 

spectroscopic interferences that are present. Other considerations include the number of dependent variables that will 
be used in the algorithm. As a general rule, the number ol samples in a learning set should be at least ten times the 
number of dependent variables employed. In the presence of known and unknown interferences, multiple sample cali- 
bration seeks to minimize the effects of interferences. Thus, the learning set solutions typify the interferences and var- 
5 iabilrty that will be experienced in actual measurements. 

Principal component analysis (PCA) is one method oi data reduction which can be performed in the application of 
chemometric techniques to spectroscopic measurement of chemical analytes in a complex matrix. PCA is used to 
reduce the dimensionality of a large number of interrelated variables while retaining the information that distinguishes 
one component from another. This reduction is effected using an eigenvector transformation of an original set ol inter- 
im related variables (e.g., an absorption spectrum) into a substantially smaller set of uncorrelated principal component 
(PC) variables that represents most of the information in the original set. The new set ol variables is ordered such that 
the first few retain most of the variation present in all of the original variables. See, e.g., Jolliffe, L.T., Principal Compo- 
nent Analysis, Sprinter- Verlag, New Vbrk (1986). More particularly, each PC is a linear combination of alt the original 
measurement variables. The first is a vector in the direction of the greatest variance of the observed variables. The suc- 
is ceeding PCs are chosen to represent the greatest variation of the measurement data and to be orthogonal to the pre- 
viously calculated PC. Therefore, the PCs are arranged in a descending order of importance. 

The term "weighting constant" is meant to include, among other techniques known in the art, the wavelength coef- 
ficients of partial least squares regression and/or principal component regression, or any constant obtained from any 
statistical calibration that can be used to calculate values (such as analyte concentration) for unknown samples. A 
20 "wavelength weighting factor" is an embedment of a weighting constant which is used in the construction of an optical 
transfer cell or optical filter means capable of extracting wavelength-specific information from spectral data. The wave- 
length-specific information can be used to determine desired values relating to the sample undergoing analysis (e.g., 
analyte concentration). A wavelength weighting factor can be embodied as a particular filter density (e.g., neutral or 
wavelength-specific), filter thickness, or the like, such parameters having been determined using the above-descrtoed 
25 statistical calibration techniques. 

The term "optical transfer ceir encompasses any optically active element that partially absorbs incident radiation 
in the visible, ultraviolet or infrared spectral regions, and wherein the absorption is selective with respect to wavelength. 
For the purposes of the present invention, an optical transfer cell generally comprises an optical filter means having 
absorption characteristics that were derived from a partial least squares or principal component regression analysis. 
30 The optical filter means is used to selectively pass one or more wavelengths from a radiation sourca said wavelengths 
having high correlation with a selected analyte concentration. "High correlation." or "close correlation" refers to the 
quantitative association between the absorption spectrum at a particular wavelength and a particular analyte concen- 
tration, wherein the two variables have a correlation coefficient (r) of 0.9 or higher. 

"Correlation spectrometry" refers to a spectrometric method wherein radiation is passed through a sample cell con- 
as taining an unknown amount of an analyte of interest, and a reference cell that contains a specified amount of the ana- 
lyte which is to be measured. By detecting the intensity of radiation which has passed though each cell, a signal which 
is representative of the difference in the radiation intensities can be generated, wherein that difference is proportional 
to the analyte concentration in the sample. 

A "negative correlation filter is provided when the sample cell used in correlation spectrometry comprises an opti- 
40 cal fflter having an absorption spectrum that is sufficient to block radiation over a selected bandwidth corresponding to 
the maximum expected concentration range of the analyte being measured. Thus, as the level of the analyte changes 
in a sample, the radiation passing through the negative correlation filter remains unaffected, as the filter has blocked all 
of the radiation in the analyte- specif ic absorption bands. 

A "positive correlation filter" is provided when the sample cell used in correlation spectrometry comprises an optical 
4S filter having an absorption spectrum sufficient to pass radiation of a particular wavelength corresponding to the target 
analyte and not to other absorbing analytes. Thus, the positive correlation filter provides an optimal transfer function 
that is highly correlated with the analyte concentration in the sample being measured. An ideal positive correlation filter 
would correlate perfectly with a target analyte (i.e., the correlation coefficient r would be +1.0), and not correlate at all 
with all other interfering absorbing analytes In a particular sample (r would be 0.0). The synthesis of positive correlation 
so filters is carried out herein using chemometric techniques to determine appropriate wavelength weighting factors as 
described above. 

A "neutral density filter" refers to a standard optical fflter means having a flat absorption spectrum. A neutral density 
filter can be used in concert with correlation filters in a filter system to provide a weighting factor to attenuate absorb- 
ances due to the analyte of interest at selected wavelengths and further improve the accuracy of the correlation pro- 
55 vided by the system. A neutral density filter can have an absorption spectrum sufficient to attenuate radiation equally at 
all wavelengths in the range of interest. A reference cell used in correlation spectrometry can comprise a neutral density 
filter. 

As used herein, an "aqueous medium" encompasses any substrate relating to. made from, or comprising water. 
Thus, an aqueous medium includes mediums wherein water is the major component, where water is present in an 
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amount of about 50%, and solutions wherein water is a solvent but is presenl in amounts of less than about 50%. Aque- 
ous mediums are specifically defined herein to include mammalian tissue. 

The term "blood analyte" refers to a blood constituent that is absorbing in the near- 1 R range, the measurement of 
which is uselul in patient monitoring or in the provision ol health care. 
5 As used herein, the term "near infrared" or "near- 1 R" encompasses radiation in a spectrum ranging from about 660 

to about 3000 nm. more preferably from about 1 050 to about 2850 nm ( and most preferably from about 1 100 to about 
2500 nm. 

The term "background absorption*' relates to the overall or base level of optical absorption of an aqueous sample 
which is to be analyzed, from which the absorption of a selected constituent departs at one or more characteristic wave- 

10 lengths to an extent indicative of the concentration of the selected constituent. When the level of background absorption 
is high in relation to the characteristic absorption of the selected constituent, such as in complex aqueous media where 
numerous Interfering constituents are found, accurate measurement of the magnitude of a slight change In the absorp- 
tion at the characteristic wavelength of a constituent of interest requires application of the chemometrics techniques 
described herein. This is particularly so in applications wherein the overall concentration of the constituent of interest is 

is low relative to the aqueous medium (e.g., in the measurement of blood analyte constituents). 

General Methods 

In complex spectral backgrounds such as those present in the spectral analysis of animal tissue or blood samples, 
20 interfering substances can produce significant measurement errors. Thus, a primary object of the present invention is 
to design an optically active element (e.g., an optical transfer cell, or positive correlation filter system) to provide the 
optimal wavelength weightings for accurate analyte estimation. Synthesis of the optimal positive correlation filter can be 
achieved according to the invention using two mathematical techniques, rotated principal components (RPC) analysis 
and least squares analysis. 

25 

A. Rotated Principal Components Analysis Method: 

RPC analysis uses an eigenvector transformation of the original spectrum to provide a set of orthogonal uncorre- 
lated principal component (PC) variables. These PCs are then individually correlated with a target analyte. Through a 
30 series of orthogonal rotations, these principal components are then transformed so that a single, rotated PC contains 
all of the information for a particular target analyte (maximum r) and the other principal component contains no informa- 
tion (r - 0). The parameter vector relating this rotate PC to the original spectrum provides the wavelength weighting fac- 
tors for the optimal positive correlation fitter. 

3S 1. Eigenvector Transformation. 

Eigenvector analysis (or principal components analysis) involves the rotation of an original data vector (X) to form 
a new data vector (X') which can be expressed as 

40 X= AX 

wherein A = an m x m matrix, and m = the number of elements in the original vector (X). 
The matrix (A) is derived from the covariance matrix or the correlation matrix of a set of representative samples 
from a data population of interest. The covariance matrix is calculated from a sample data set. This covariance matrix 
45 is then transformed by a process of diagonalization into a new PCs matrix in which the rows of the matrix are ordered 
by a set of roots known as eigenvalues. The horizontal vectors of this new matrix (A) are called eigenvectors. The sub- 
ject matrix can be used to transform an original spectral data vector into a transformed vector of independent uncorre- 
tated variables. In this regard, the concept of independence is critical in near-IR spectral analysis of a complex chemical 
sample (such as blood) In order to remove the effects of Interfering chemical analytes. Uncorrected variables that relate 
so to particular analytes serve as the vehicles of accurate estimations of the analyte concentration. 

The above-described transformation does not provide uncorrected variables that are specific to a particular ana- 
lyte. Rather, it generates uncorrelated variables of decreasing levels of variance that correlate to varying degrees with 
a number of different analytes. 

55 2. Principal Component Rotation for Maximal An alyte Correlation. 

Typically, two to four of the PCs derived from the above PC analysis account for 85 - 98% of the variance in the orig- 
inal variables. These high variance variables are not necessarily strongly correlated with a particular analyte. Other 
analytes or matrix effects may contrtoute greater changes in absorbances. A correlation matrix of a dependent variable 
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analyte provides a convenient basis to rotate a series of PC pairs to produce a single revised PC thai more completely 
"explains" a particular analyte. The revised PCs that relate to a particular analyte will thus supply the parameters for 
accurate estimation ol analyte concentrations. 

Each pair of PCs defines a plane. The two PCs that define the plane can be rotated with the purpose f converting 
s the two PCs into a single PC with all of the information content of the original PC pair. One method of completing this 
conversion is to implement the following computational procedure. 

First, the two highest correlated PCs tor a particular analyte are selected. This PC pair i6 rotated until one PC con- 
tains all of the information content of the original pair. The rotation angle may be iteratively determined using an exhaus- 
tive experimental sequence or it may be analytically determined without iteration. Thus, 

w 2 2 

for PCYR', +R 2 .and 



forPC' a : R 2 =0. 

is Continuing with the procedure, the next highest correlated PC (PC3) is selected and rotated with PO\ until, 

fbrPCgiR'g 2 - R 3 2 + R , 1 2 .and 



for PC 1 .,: R't =0. 



20 



The above rotational sequence is repeated until all ol the analyte correlating vectors have been treated. The single, 
rotated PC generated by the rotational sequence contains all of the information for a particular target analyte (maximum 
r) and the other PC contains no information (r = 0). Further, the subject PC is independent of the absorbance contribu- 
tions produced by other optically active analytes in the sample. The.parameter vector relating this rotated PC to the orig- 
25 inal spectrum provides the wavelength weighting factor for the optimal positive correlation filter (optical transfer cell). 

B. Least Squares Analysis Method : 

In the application of the least squares method, the optimal transfer lunction is obtained directly. Given a collection 
30 of vectors. 

X|, wherein / ■ 1, n, 

a vector (v) must be found such that x^v wilt lie as close as possible to a line when x?v is plotted against c, (the con- 
35 centration value for x,). Thus, a and b must also be found, such that 

ax, T v + b = c j . for all /'. 

Since equality cannot be actually attained, a minimum must be accepted, in a least squares sense, of 

40 

|1axj v + b-Cj|,for all /. 
Thus, the problem is formulated as follows, A matrix (A) is formed: 



A- 



50 



X \ 



\x n T ) 



ss 'and vectors 
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c = 




b = 






• 







w When the objective function is rewritten as 

mln (a. 6. v) = ||aAv + b • c||, and 
by letting v = av , the above function decreases to 

75 

min v, b ||Av + b - c||. 

Now, defining 



25 the problem becomes 

min(v) = ||Av - c||. 
This relationship is well known and has the solution 

30 

V = (A T A)* 1 A T C 

The vector (v) is thus the optimal weighting function to implement the optical transfer function in the positive correlation 
filter. 

35 Determined over a large number of samples, both the rotated principal components and the least squares analysis 
methods will likely produce substantially similar optical transfer functions. However, for more limited sample sets, the 
rotated principal components methods will be more effective due to the orthogonal nature of the optimal transfer func- 
tion. The transfer functions are then used to determine the parameters for the positive correlation filters of the present 
invention. 

40 The aforementioned positive correlation filters are employed herein in a variety of devices configured to extract ana- 
lyte concentration information from complex spectral media using near-IR spectrometer devices. A number of near-IR 
spectrometers, Fourier Transform near-IR spectrometers, and modified 1R spectrometers of conventional design can be 
used in the practice of the invention. The preferred mode of operation is reflectance. Suitable spectrometers include the 
NIR-Systems model 5000 analyzer, available from Perstorp Analytical (Silver Spring, MD), and the Luminar 2000 ana- 

45 lyzer, available from Brimrose (Baltimore, MD). The techniques for sampling, measuring, and signal processing can be 
conventional and are well known to those of skill in the art. 

Referring now to Figure 1 . an apparatus for measurement of analyte concentration in a complex aqueous medium 
is generally indicated at 1 0. The apparatus includes a radiation source 1 2 which is capable of generating a spectrum of 
near Infrared (near-IR) radiation within the range of approximately 600 to approximately 2500 nm. A number of suitable 

so radiation Sources are known in the art. such as incandescent light sources directed across interlerence filters, halogen 
light sources (modulated by an associated chopper wheel), laser light sources, laserdiode arrays, or high speed light- 
emitting diode (LED) arrays. 

Radiation from the source 1 2 is directed to sample interface optic means 14 which reflectively transmits the radia- 
tion to a medium 16 containing the Wood analyte of interest. After the radiation has interrogated the medium, spectrally 
55 modified radiation is transmitted back to the apparatus 10 through further sample interface optic means 18. The appa- 
ratus 10 can be designed to closely interface with the medium 1 6. such as where the radiation is directed to the medium 
and modified radiation transmitted back to the apparatus using an optically active means 20 (e.g.. light converging 
means such as a lens, or beam deflection optics) known in the art. Alternatively, the sample interface optic means 14 
and 16 can comprise ffoer optic waveguides coupled to the apparatus in order to enable more remote medium interfac- 
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ing. Other configurations are provided wherein a single fiber optic bundle is employed to transmit the radiation to and 
from the medium. An optrode disposed at the end of the single bundle transmits the near-IR radiation to the sample 
medium 16 and receives spectrally modified radiation therefrom which is directed back through the bundle to the appa- 
ratus 10. Sapphire or high-grade quartz can be used as optical elements in the above fiber optic waveguides, as those 

5 materials have very good transmission characteristics in the near-IR spectral range. Other suitable fiber optic 
waveguides include fluoride-based systems known in the art. 

In each of the above-described configurations, radiation from sample interface optic means 14 is directed to the 
medium 16 where the near-IR radiation interacts with the constituents of the sample and is absorbed, dispersed, dif- 
fused and reflected. This spectrally modified radiation is then collected and directed to the apparatus 10. 

jo The spectrally modified radiation transmitted by sample interface optic means 1 8 is then directed to a means 22 for 
receiving and directing the radiation into a beam path, such as a coilimating lens or the like. The beam path communi- 
cates with means 24 for passing selected wavelengths, such as a diffraction grating system or selective filtering means. 
Use of a diffraction grating system serves to eliminate wavelengths not having a relationship to the analyte of interest, 
and divides the radiation into distinct wavelengths. Radiation passing through the means 24 is directed through a beam 

is splitter 26, whereby the radiation is divided into two beams, respectively indicated at 28 and 30. Suitable beam splitter 
means include mirror configurations such as those known in the art. 

A first beam 28 from beam splitter 26 is communicated with an optical transfer cell 32 which is disposed in the beam 
path. The optical transfer cell 32 comprises a positive correlation filter system having absorption spectra sufficient to 
accept the spectrally modified near-IR radiation from the beam 28 and selectively pass one or more wavelengths there- 
to from having high correlation with a selected analyte concentration. The positive correlation filter system includes one or 
more optically active filter means having absorption characteristics determined using a partial least squares or principal 
component regression of an original spectrum of a sample containing an analyte of interest. 

Positive correlation filter systems can be constructed using a suitable substrate layer that is capable of transmitting 
radiation at least in the 1 100-2500 nm range. The substrate layer is generally ooated with one or more layers of metals 

25 and/or oxides that are conventional in the art. Such coatings can be applied to the substrate using emulsion or chemical 
vapor deposition (CVD) techniques well known in the art. In preferred embodiments, the positive correlation filter sys- 
tem comprises a photographic mask with spectral lines of optical density proportional to an optimal weighting function 
determined using rotated principal components or least squares analysis methods. One particular filter system has a 
filter density and/or thickness that comprises a wavelength weighting factor means, wherein the weighting factor pro- 

30 vides enhanced positive correlation of a passed wavelength with the concentration of an analyte in the selected sample 
medium. The weighting factor means is an embodiment of a weighting constant determined using chemometrics, 
wherein the constant allows for the construction of an optical filter that is able to extract wavelength-specHic information 
from spectraJ data. The wavel ength -specif ic information can then be used to readily determine analyte concentration. 
Thus, the optical transfer cell 32 allows for the measurement of analyte concentration in a variety of complex medi- 

35 urns, such as any aqueous medium having a complex spectral background (e,g., in which interfering constituents in the 
medium can produce significant measurement errors). Particularly, the optical transfer cell 32 finds use in applications 
where there exists significant background absorption from which the absorption of a specific constituent of interest 
departs at one or more characteristic wavelengths to an extent indicative of the concentration ot that constituent. When 
the level of the constituent of interest is low relative to the medium, the need for such an optical transfer cell becomes 

40 more acute. 

In one application, the optical transfer cell 32 can be used in the determination of Wood analyte concentrations. The 
blood analyte can be present in an in vitro sample medium (e.g. ( a blood sample), or can be present in an in vivo tissue 
sample. In particular applications, the Wood analyte of interest is a Wood electrolyte (e.g., calcium, potassium, sodium, 
chloride and bicarbonate (COJ), blood hydrogen ion concentration (pH). an organic compound present in blood such 

45 as glucose and urea (BUN), or a Wood constituent such as hematocrit or hemogloWn. In each of the aforementioned 
devices, the optical transfer cell 32 includes a positive correlation filter system that has been synthesized tor the deter- 
mination of the concentration of the specific analyte of interest. As will be readily appreciated by those of skill in the art, 
the optical transfer cell can be detachabry mounted in the apparatus 10, whereby a variety of different cells can be cas- 
setted into and out of the apparatus when one is desirous of quantifying different analytes. 

so Referring still to Figure 1 . the correlated wavelengths passed by optical transfer cell 32 are received by a detection 
means 34 for converting the optical information (correlated wavelengths) into a signal that is representative of the inten- 
sity of that information. In one particular apparatus, the detection means 34 comprises a lead sulfide (PbS) photodetec- 
tor able to scan the range of wavelengths from about 1 1 00 to at least about 2500 nm in steps of 1 nm. Such detection 
means convert the passed radiation (wavelengths) into signals representative of the intensities of those wavelengths. 

55 A second beam 30 from beam splitter 26 is communicated to a further detection means 36, whereby the spectral 
information is converted into signals representative of the intensity thereof. The second beam 30 of spectrally modified 
radiation may pass through a reference optical transfer cell 38 disposed in the path of said beam. The reference cell 38 
can comprise a neutral density fitter means having absorption characteristics sufficient to attenuate radiation equally 
over a selected range of near-IR wavelengths. The signals from detection means 34 and 36 are communicated to a 
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means 40 for converting those signals into digital signals, e.g. digital signals indicative of the differences between the 
correlated wavelengths emerging from the positive correlation filter and the attenuated wavelengths emerging from the 
neutral density filter. More particularly, the signals are communicated to an analog/digital (A/D) converter for conversion 
to digital form. The digital information is then readily input into a microprocessor 42 to calculate actual analyte ooncen- 

5 tration. The concentration can conveniently be visualized through a digital display means 44. 

Various optional components for use in the apparatus 10, such as first and second radiation focusing means, 
respectively indicated at 46 and 48 and capable ol receiving and directing the radiation emerging from optical transfer 
cells 32 and 38 into a beam path, will be apparent to those skilled in the art. In some apparatuses, such focusing means 
can comprise a focusing lens or the like. Alternatively, the apparatus 10 can include, for example, a single detection 

10 means capable of receiving spectral information from the split beams 28 and 30. 

Referring now to Figure 2. an alternative apparatus for measurement ol analyte concentration in a complex aque- 
ous medium is generally indicated at 60. The apparatus includes a radiation source 62 which Is capable ol generating 
a spectrum of near infrared (near-IB) radiation within the range of approximately 660 to approximately 2500 nm. Radi- 
ation from source 62 is transmitted to an optically active means 64 for receiving and directing the radiation into a beam 

is path and/or for passing selected wavelengths, such as a collimating lens, selective filtering means or the like. In one 
case, the optically active means 64 can comprise a monochromator, wherein a control signal obtained from a micro- 
computer is used to step the position of the monochromator relative to the beam path. 

Near-IR radiation emerging from means 64 is communicated through a beam splitter 66. whereby the radiation is 
divided into two beams, respectively indicated at 68 and 70. The first beam 68 from beam splitter 66 is transmitted to a 

20 sample medium 72 containing an unknown concentration of an analyte of interest. In Figure 2. the sample medium 72 
comprises a sample cell formed from a suitable substrate capable of transmitting radiation in the near-IR range of inter- 
est. In one case, the sample can comprise a Wood serum sample, wherein it is desired to determine the concentration 
of Wood anatytes. Alternatively, the first beam 68 can be transmitted to a sample surface such as a tissue surface using 
direct interfacing means or indirect interfacing means (e.g., fiber optic waveguides means) as have been described 

25 supra. In this manner, the concentration of blood analytes present in the tissue sample can be non-invasively deter- 
mined using reflective neaMR measurement of the absorption spectrum of racfiation which has interacted with the tis- 
sue sample. 

Spectrally modified radiation that has interacted with the constituents of the sample (including the analyte of inter- 
est) is then collected and directed to an optical transfer cell 74 which is disposed in the beam path. The optical transfer 
30 cell 74 comprises a positive correlation filter system having an absorption spectrum sufficient to accept the radiation 
and selectively pass one or more wavelengths therefrom having high correlation with the concentration of an analyte of 
interest, and substantially no correlation with interfering components present in the sample. Radiation emerging from 
the optical transfer cell 74 is received by a detection means 76 for converting the spectrally modified radiation into a sig- 
nal that is representative of the intensity of that radiation. The detection means can comprise a broad spectrum photo- 
ns detector such as an PbS photodetector or the like. 

Referring still to Figure 2, the second beam 70 from beam splitter 66 is transmitted to an optically active element 
78 which is disposed in the beam path. In one configuration, the optically active element 78 comprises a neutral density 
filter means having absorption characteristics sufficient to attenuate radiation equally over a selected range of near-IR 
wavelengths. In an alternative configuration, the optically active element 78 is an optical transfer cell that includes a 
40 positive correlation filter system having an absorption spectrum that is identical to the absorption spectrum of optical 
transfer cell element 74. The radiation emerging from the optically active element 78 is received by a detection means 
80 for converting the radiation into a signal that is representative of the intensity thereof. 

The positive correlation filter systems can be constructed as described above. Thus, the subject filters can com- 
prise a single substrate layer having an optically active coating which imparts absorption characteristics capable of 
45 selectively passing one or more wavelengths having high correlation with a particular analyte concentration. In pre- 
ferr d embodiments, the positive correlation filter system comprises a plurality of filter layers, each layer having a 
selected filter density and/or filter thickness suitable to provide a desired absorption characteristic. In one case, at least 
one layer of the system has a filter density and/or thickness that comprises a wavelength weighting factor means, 
wherein the weighting factor provides enhanced positive correlation of a passed wavelength with the concentration of 
so an analyte in the selected sample medium. 

The signals generated by detection means 76 and 80 are then communicated to a means 82 for converting those 
signals into a digital signal indicative of the ratio of the intensity of the radiation emerging from the source 62 and the 
corresponding spectrally modified radiation emerging from the sample. In this manner, changes in the intensity of the 
radiation emerging from the source 62 can be corrected for. eliminating a potential source of error in the measurement 
55 obtained from the system. Further, the ratio of the signals can then be converted to digital form and interpreted to deter- 
mine analyte concentration using an internal microprocessor 84 system or an associated system using methods known 
in the art. 

In particularly preferred embodiments, the microprocessor calculates analyte concentration by application of a che- 
mometrics algorithm to the ratio signal. The appropriate algorithm can be determined using the above-described che- 
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mometrics techniques, such as a least squares analysis or totaled principal component analysis of an original 
absorption spectrum of the analyte ol interest 

In a further aspect of the invention, a non-invasive method of measuring the concentration of a blood analyte is pro- 
vided. The method generally employs a near-IR spectrometer apparatus such as those described above. Particularly, 

5 a method is provided wherein near-IR radiation is directed to a medium containing a Wood analyte of interest. In one 
case, the radiation is caused to impinge upon the skin surface of tissue, such as upon a subject's forearm. The radiation 
is directed at an angle toward the tissue such that it is absorbed by the tissue material near the surface and reflected 
as diffuse radiation. In this manner, the radiation is spectrally modified as a result of infrared absorptions by the blood 
and tissue constituents. Portions ol the interrogating near-IR radiation are absorbed, dispersed, diffused and reflected 

to from the Wood constituents contained within the tissue source. This spectrally modified radiation contains information 
specific for each optically active Wood constituent. However, due to the complexity of the chemical and physical matrix 
of the tissue, a large amount of background spectral Information is also obtained which must be removed or accounted 
for in order to extract analyte-specific information. The near-!R radiation used to interrogate the medium preferably 
includes radiation in the range of approximately 1 1 00 to 2500 nm. 

75 In a further preferred aspect of the invention, a non-invasive method is provided for determining blood glucose lev- 
els using the above-described transcutaneous forearm measurement techniques. The vibrational motion of the blood 
glucose molecules can be measured using diffuse- reflective near-IR radiation, wherein the vibrational motion includes 
both rotational and translations! motion, including overtone vibrations and combined vibrations. Of these vibrations, the 
overtone vibrations are dominant. 

20 Such transcutaneous, reflective measurements carried out in the forearm area provide a more uniform representa- 
tion of Wood glucose than prior methods that are carried using transmissive near-IR finger measurements. Near-IR light 
transmission in the finger area interacts with a mixture of arterial, venous and capillary blood, as well as bone, cartilage 
and other subcutaneous material. On the other hand, reflective arm measurements are limited to capillary Wood (and 
interstitial fluid) almost exclusively. Since capillary, venous and arterial Wood glucose concentrations are known to differ 

25 slightly, the measurement of a single blood state (e.g., in transcutaneous forearm measurements) provides lor a more 
consistent monitor of blood glucose. Reflective arm measurements are also much less sensitive to sensor location 
placement than transmissive finger measurements, wherein variation of sensor placement and rotation can lead to 
errant spectral measurements. Furthermore, skin pigmentation variations from racial and individual differences cause 
absorption variations in the transmissive near-IR region (e.g., 800 to 1000 nm), but not in the reflective near-IR region 

30 (e.g., over 1000 nm). 

In each of the above methods, the near-IR radiation that was directed to the medium containing the analyte ot inter- 
est is caused to be reflected by the medium to provide a spectrally modified beam of radiation. This spectrally modified 
radiation is directed into a path through an optical transfer cell disposed within the beam path. The optical transfer cell 
comprises a positive correlation filter system, as has been described in detail above. 

35 Accordingly, the invention provides a positive correlation fflter system configured to accept the beam of spectrally 
modified near-IR radiation and selectively pass one or more wavelengths therefrom. The passed wavelengths are 
selected to have high correlation with an analyte concentration. The methods further entail measuring the intensity of 
the passed wavelengths and generating a signal representative of the intensity of the wavelengths. The resultant signal ' 
can then be converted into an indication of analyte concentration using methods known in the art. In a preferred emtood- 

40 iment. the analyte concentration is determined using an application of a chemometrics algorithm to the resultant signal. 
It is to be understood that while the invention has been described in conjunction with preferred specific embodi- 
ments thereof, the foregoing description, as well as the examples which follow, are intended to illustrate and not limit 
the scope of the invention. Other aspects, advantages and modifications within the scope of the invention will be appar- 
ent to those skilled in the art to which the invention pertains. 

45 

Example 1 

Transcutaneous near-IR reflective measurement of blood analyte concentrations was carried out in parallel with 
actual laboratory testing of patient Wood samples to evaluate the accuracy of such methods and to develop training set 

so data for the above described chemometric techniques. More particularly, a near-IR photodiode array analyzer operating 
reflectively in the 1 100-1800 nm range was used to test two separate patient populations (at the Bethesda Naval Hos- 
pital (Maryland) and the Froedtert Lutheran Memorial Hospital (Milwaukee), respectively). The testing apparatus 
employed a halogen light source that was coupled to an optical fiber means through a shutter and attenuator means. 
The shutter means allowed for light Wocking sufficient to provide for a dark reading of the associated photodetector 

55 array. Near-IR radiation in the 1 100-1800 nm spectral region was transmitted to a reflective optical probe that was cou- 
pled to the patient sWn surface (in the forearm area). Diffuse back- scattered radiation that had penetrated the epidermis 
and dermis was collected by a receiving optical fiber means that transmitted the radiation to a spectrograph where it 
was detected by an indium gallium arsenide photodeiector array (256 x I) and converted to digital form by array inter- 
face electronics. The spectral signals were then interfaced to a microcomputer for signal processing and pattern recog- 
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nition. 

The results of the non-invasive transcutaneous analysis of the following blood electrolyte analytes: calcium (CA*). 
potassium (K*) t sodium (NA*), chloride (CL') and bicarbonate (C0 2 ); Wood glucose (GLU), Wood urea (BUN); hemat- 
ocrit (HCT) and hemoglobin (HEM) are set forth in Tab! I. An effective measurement is characterized by a small aver- 
5 age error, a slope (b) approaching 1 .0 and a high T- value (T). The slope is defined as the tangent of the least squar fit 
line of predicted values/actual values. The predicted values were obtained using the transcutaneous measuring appa- 
ratus of the invention and the actual values were obtained by physical laboratory testing of patient blood samples. The 
T-value can be described as a "tracking value" that indicates how well the predicted value tracked with changes in the 
actual value. The T-value is defined as the slope/standard error of the slope. 
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TABLE I 



75 


Analyte 


Mean Value 


Range 


Average 
Error 


(b) 


T 


CA + mg/dl 


9.63 mg/dl 


8.5-10.6 


0.21 mg/dl 
2.18% 


0.57 


7.01 


20 


K + mmol/L 


4.44 mmol/L 


3.5-5.5 


0.16 mmol/L 
3.60% 


0.75 


10.82 




NA + mmol/L 


140.51 mmol/L 


135-148 


1.24 mmol/L 
0.88% 


0.59 


6.95 


25 


CL" mmol/L 


103.68 mmol/L 


94-110 


1.13mmol/L 
1.09% 


0.57 


9.56 




CO2* mmol/L 


28.56 mmol/L 


22-35 


1.55 mmol/L 
5.43% 


0.46 


4.12 


30 


QLUmg/dl 


98.07 mg/dl 


60-115 


13.97 mg/til 
14.24% 


1.00 


14.61 




BUN mg/di 


15.58 mg/dl 


7-26 


1.99 mg/dl 


0.67 


9.17 


3S 








12.77% 






HCT% 


41.78% 


33-55 


2.62% 

8.27% of mean 


0.41 


5.94 




HEM g/dl 


14.32 g/dl 


12-18 


0.89 g/dl 


0.48 


6.45 


40 








6.22% 







45 Example 2 

The Transcutaneous near-IR reflective measurements of blood analyte concentrations carried out in Example 1 
were repeated using a additional 500 patients at the Naval Hospital in San Diego, CA. The data generated were used 
10 carry out a training set/test set sequence using the chemomeulc methods of the present Invention. Trie results of a 

50 "leave-one-out" principal component regression analysis of the following blood electrolyte analytes: calcium (CA*), 
potassium (K + ), sodium (NA+), chloride (CL ) and bicarbonate (C0 2 ); Wood glucose (GLU). blood urea (BUN); hemat- 
ocrit (HCT) and hemoglobin (HEM) are set forth in Table II. Once again, an effective measurement is characterized by 
a small average error, a slope (b) approaching 1 .0 and a high T-value (T). The predicted values were obtained using the 
transcutaneous measuring apparatus of the invention and the actual values were obtained by physical laboratory test- 

55 ing of patient blood samples. 
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Claims 



An optical transfer cell for use in the measurement of analyte concentration in a complex aqueous medium, com- 
prising a positive correlation filter means adapted to accept near infrared radiation and to selectively pass a wave- 
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length from said radiation having high correlation with the analyte concentration. 

2. The optical transfer cell of claim 1 , wherein the medium comprises body tissue and the analyte comprises a Wood 
analyte. 

5 

3. The optical transfer cell of claim 2, wherein the Wood analyte is selected from the group consisting of calcium, 
potassium, sodium, chloride, bicarbonate (COJ, hydrogen ion concentration (pH). glucose, urea (BUN), hematocrit 
and hemoglobin. 

10 4. The optical transter cell of claim 3. wherein the near infrared radiation includes wavelengths in the range of approx- 
imately n 00 to 2500 nm. 

5. The optical transfer cell of claim 1 , wherein the positive correlation filter means is adapted to selectively pass a plu- 
rality of wavelengths, each having a high correlation with the analyte concentration. 

75 

6. The optical transfer cell of claim 1 , wherein the positive correlation filter means further comprises a plurality of lay- 
ers, each layer having a selected near infrared absorption characteristic such that said fitter means passes a pop- 
ulation of wavelengths, each having a high correlation with the analyte concentration. 

20 7. The optical transter cell ol claim 6, wherein at least one layer further comprises wavelength weighting factor means, 
said weighting factor providing enhanced positive correlation of the passed wavelength to the analyte concentration 
in a selected medium. 

8. The optical transfer cell of claim 7, wherein said weighting factor is determined using rotated principal component 
25 analysis of an original absorption spectrum of the analyte. 

9. The optical transter cell of claim 7, wherein said weighting factor is determined using a least squares analysis of an 
original absorption spectrum of the analyte. 

30 10. An apparatus for measurement of analyte concentration in a complex aqueous medium, said apparatus compris- 
ing: 

a source of near infrared radiation; 

means for reflectively transmitting the radiation to the aqueous mecSum and retransmitting spectrally modified 
as radiation reflected from said aqueous medium; 

means for receiving and directing the spectrally modified radiation into a beam path; 
an optical transfer cell disposed in the beam path, said cell comprising positive correlation filter means adapted 
to accept the radiation beam and selectively pass a wavelength from said beam having high correlation with 
the analyte concentration; 

40 means for receiving the passed wavelength and for converting the same into a signal representative of the 

intensity of said wavelength; and 
means for calculating analyte concentration from the signal. 

11. The apparatus of claim 1 0. wherein the means for calculating analyte concentration comprises application of a che- 
45 mometrics algorithm to the signal. 

1 2. The apparatus of claim 10. wherein the positive correlation titer means further comprises a plurality of layers, each 
layer having a selected absorption characteristic whereby said filter means passes a population ol wavelengths 
having high correlation with the analyte concentration. 

50 

13. The apparatus ol daim 12, wherein at least one layer of the filter means further comprises wavelength weighting 
factor means, said weighting factor providing enhanced positive correlation of the passed wavelength to the analyte 
concentration in a selected medium. 

55 14. The apparatus of claim 13, wherein the medium comprises body tissue and the analyte comprises a blood analyte. 

15. The apparatus of claim 1 4, wherein the blood analyte is selected from the group consisting of calcium, potassium, 
sodium, chloride, bicarbonate (CO2). hydrogen ion concentration (pH), glucose, urea (BUN), hematocrit and hemo- 
globin. 
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16. The apparatus of claim 14, wherein the near infrared radiation includes wavelengths in the range of approximately 
1100 to 2500 nm. 

17. A non-invasive method of measuring the concentration of a blood analyte. comprising: 

5 

(a) directing near infrared radiation to a medium containing the blood analyte; 

(b) causing the radiation to be reflected by the medium to provide a spectrally modified beam; 

(c) directing the modified beam into a path through an optical transfer cell disposed in the beam path, said cell 
comprising positive correlation filter means adapted to accept the beam and selectively pass a wavelength 

to from said beam having high correlation with the analyte concentration; 

(d) measuring the intensity of the passed wavelength and generating a signal representative of the intensity of 
the wavelength; and 

(e) converting the signal into an indication of analyte concentration. 

is 18. The method of claim 1 7, wherein the medium comprises an aqueous solution. 

19. The method of claim 1 7, wherein the medium comprises a solid substrate. 

20. A non-invasive method of measuring the concentration of a blood analyte, comprising: 

20 

(a) directing near infrared radiation to the surface ol a substrate containing the blood analyte; 

(b) causing the radiation to be reflected by the surface to provide a spectrally modified beam; 

(c) directing the modified beam into a path through an optical transfer cell disposed in the beam path, said cell 
comprising positive correlation filter means adapted to accept the beam and selectively pass a wavelength 

25 from said beam having high correlation with the analyte concentration; 

(d) measuring the intensity of the passed wavelength and generating a signal representative of the intensity of 
the wavelength; and 

(e) converting the signal into an indication of analyte concentration. 

30 21 . The method of claim 20. wherein the substrate comprises human skin or mucosal tissue. 

22. The method of claim 1 7, wherein the step of converting the signal into an indication of analyte concentration com- 
prises application of a chemometrics algorithm to the signal. 

35 23. The method of claim 17, wherein the blood analyte is selected from the group consisting of calcium, potassium, 
sodium, chloride, bicarbonate (C0 2 ), hydrogen ion concentration (pH), glucose, urea (BUN), hematocrit and hemo- 
globin. 

24. The method of claim 1 7. wherein the near infrared radiation includes wavelengths in the range of approximately 
40 1100 to 2500 nm. 
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